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MF  Radio  Field  Strength  Measurements 
in  Desert  Terrain  Near  Yuma,  Arizona 


I.  INTRODUCTION 

The  possibility  that  a  network  of  underground  tunnels  or  silos  containing  the 
MX  mobile  missile  might  be  constructed  in  the  desert  regions  of  the  southwestern 
United  States  has  renewed  interest  in  techniques  for  survivable  communications 
between  buried  terminals  which  may  be  up  to  50  km  apart  in  adjacent  valleys  sepa¬ 
rated  bv  high  mountain  ridges.  One  proposed  method  is  to  use  an  MF  radio  link 
with  the  transmitting  and  receiving  antennas  located  on  the  underground  mobile 
launch  control  centers.  The  propagation  path  for  such  a  link  would  include  two  seg¬ 
ments  through  the  trench  wall  and  soil  overburden  in  addition  to  that  over  the  desert 

terrain.  Although  MF  propagation  curves  over  homogeneous  earth  are  readily 
1  -4 

available,  propagation  anomalies  caused  by  topographical  and  electrical  irregu¬ 
larities  unique  to  the  proposed  MX  sites  need  to  be  studied  in  situ. 


(Received  for  publication  24  November  1980) 
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2.  Wait,  J.  R.  ,  and  Campbell,  L.  L.  (  1953)  Transmission  curves  for  ground  wave 

propagation  at  low  radio  frequencies,  Radio  Physics  Laboratory  Report  R-l, 
Defence  Research  Telecommunications  Establishment,  Ottawa. 
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2.  EXPERIMENTAL  PROCEDURE 

Between  11  September  and  24  September  1977  personnel  from  Rome  Air 
Development  Center  (RADC/FEP)  made  field  strength  measurements  of  the  Yuma, 
Arizona  broadcast  station  KULU  (560  kHz,  1  kW,  omni  -directional  in  the  daytime) 
at  several  hundred  positions  on  Luke  A  F  Bombing  and  Gunnery  Range,  a  desert 
region  southeast  of  Yuma.  At  each  site  the  surface  magnetic  field  was  measured 
using  a  Radio  Interference-Field  Intensity  Measuring  Equipment,  Singer  Model 
NM-25T,  with  a  shielded  15-in.  -diameter  loop  antenna  oriented  for  maximum 
signal.  At  some  of  the  sites  the  loop  also  was  lowered  2  m  into  a  hole  drilled  by 
a  power  auger  to  probe  the  radio  wave  attenuation  with  depth.  The  measurements 
were  conducted  only  during  daylight  hours  to  eliminate  the  necessity  of  correcting 
for  KBLU  day-to-night  antenna  pattern  changes  and  to  minimize  skv  wave  inter¬ 
ference. 

The  measurement  sites  were  accurately  located  in  almost  every  case  by 
choosing  positions  near  the  numerous  L'SGS  markers  found  alongside  Jeep  trails. 

The  sites  are  shown  in  Figure  1,  along  with  the  location  of  the  KBLU  transmitting 
antenna  in  downtown  Yuma.  The  distribution  of  the  markers  was  such  that  no  series 
of  measurements  was  strictly  aligned  along  a  radial  path  from  KBLU.  However, 
six  sectors  were  chosen  within  which  the  terrain  could  be  characterized  reasonably 
well  by  a  single  profile,  thus  providing  a  convenient  wav  to  group  and  display  the 
data.  Dashed  lines  are  used  in  Figure  !  to  show  the  sector  boundaries.  A  repre¬ 
sentative  path  profile,  corrected  for  the  standard  tropospheric  refractive  index 
gradient3  (4/2  earth  radius),  was  constructed  for  each  sector  using  USCiS  maps. 

3.  FLAT  EARTH  PROPAGATION  THEORY 

Over  a  flat,  homogeneous,  and  well-conducting  earth  the  vertical  electric  field 
strength-distance  product  is  giver,  by 

'  E  ■  x  I  -  300  •  |F<p(i)|  (Volts)  (1) 

where  P  is  the  effective  radiated  power  in  kW, 

F(pp)  =  1  -  i  -  P  C  •  erfe(i  ^p  ) 

n.  Reference  Patti  for  Radio  Engineers,  5th  ed.  ,  (1957)  ITT  Corp.  ,  New  York 
p  74  1  ff.  B 
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F(p  )  is  the  Sommerfeld  propagation  factor,  p  is  the  numerical  distance, 
e  - 1 2  e 

t  =  8.  854  X  10  Vim,  <x  is  the  radian  frequency,  a  is  the  surface  conductivity 

(Siemens /meter)  and  kx  is  the  distance  in  wavelengths.  Eq.  (1)  can  also  be  expressed 

log  |  E  -  x  I  =  log  If!  +  log  300  ,  (2) 

where  log  =  logjg.  |  F  |  is  plotted  us  a  function  of  pg  on  a  logarithmic  scale  in 
Figure  2.  For  numerical  distances  of  the  order  of  unity  or  less, 

log  !  f!  :  -0.  183  t  pe  1  .  (3) 

At  560  kHz,  Fq.  (3)  is  valid  for  distances  out  to  at  least  50  km  if  o  =  0.  01  S/m  and 
to  100  km  foro  =  0.02  S/m.  Combining  Eqs.  (2)  and  (3), 

0.  09  15  £  o)k 

log  1  E  -  x!  -  -  - g— ^ -  x  +  log  300  .  (4) 

Equation  (4)  is  in  the  standard  form  y  =  mx  +  b,  where,  upon  the  application  of 
linear  regression  analysis  to  the  measured  (E*  x|  data,  the  slope 

i  uj  k 

111  =  -0.  0915  —  (5) 


yields  an  estimate  foro,  and  the  v-intercept 

b  =  log  300 \V  (6) 

indicates  an  apparent  radiated  power 


F'igure  2.  Amplitude  of  the 
Propagation  Factor 
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Figure  5.  Wellton  Hills 
(107.5°  -  110°)  Sector 
Path  Profile  (lower  portion) 
and  Signal  Strength-Distance 
Product 


Figure  6.  Coyote  Peak 
(110°  -  116°)  Sector  Path 
Profile  (lower  portion)  and 
Signal  Strength-Distance 
Product 
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4.  RESULTS  OF  REGRESSION  ANALYSIS 


4.1  Estimate  of  Conductivity 

The  upper  portions  of  Figures  3  to  8  show  measured  values  of  !  F  •  x |  on  a 
logarithmic  scale  vs  distance  from  the  KBLU  antenna.  The  lower  portions  are 
corresponding  terrain  profiles  which  reveal  long  flat  path  segments  occasionally 
interrupted  by  relatively  steep  ridges.  The  effective  conductivity  of  these  flat  seg¬ 
ments  was  estimated  by  fitting  a  regression  line  to  selected  data  points  not  too  close 
to  the  mountain  ridges,  and  then  using  Eq.  (5)  to  estimate  o.  The  regression  line 
is  shown  in  the  figures  as  a  solid  line,  and  the  o  estimates  are  given  in  Table  1. 
Assuming  the  data  is  normally  distributed  about  the  regression  line,  we  can  assign 
confidence  limits  to  the  calculated  values  of  u.  The  results  of  choosing  a  115  percent 
interval  are  shown  in  the  last  column  of  Table  1.  Evidently  only  the  100°  -  107.  5° 
sector  contained  enough  data  over  a  sufficiently  long  and  flat  range  to  permit  a 
reasonably  accurate  estimate  for  o,  which  was  0.022  S/m.  This  value  is  somewhat 
higher  than  the  0.  01  S/m  estimated  by  Morgan1  or  the  0.  008  S/m  estimated  bv  the 
FC'C.  8 

Some  of  the  data  scatter  within  the  flat  segments  was  due  to  the  fact  that  the 
measurements  were  not  along  radial  paths,  and  each  data  point  within  a  given 
sector  actually  represents  propagation  over  a  slightly  different  path.  Additional 
perturbations  were  produced  by  local  anomalies  such  as  ore  deposits,  faults,  power 
and  telephone  lines,  and  railroads,  as  well  as  by  focussing,  diffraction,  and  multi- 
path  effects. 


Table  1.  F.stimated  Effective  Conductivity 


Sector 
( Degrees  | 

No.  Data 
Points 

.. 

Correlation 

Coefficient 

Est  (7 
[S/m  j 

05%  Coni'.  Limits 
[  S  /  m  ) 

90  -  100 

24 

0.  39 

0.  023 

* 

100  -  107.5 

44 

0.  70 

0.  022 

0.  016  -  0.  033 

107. 5  -  110 

19 

0.  42 

0.  029 

* 

110  -  116 

21 

0.  22 

0.  026 

& 

116  -  120 

6 

0.  22 

o 

o 

o 

* 

123  -  137 

20 

0.  52 

0.  023 

0.  011  -  0.  122 

The  data  scatter  was  too  great  to  permit  an  accurate  determination  of  a. 


(Duo  to  the  number  of  references  cited  above,  they  will  not  be  listed  here,  bee 
References,  page  20.  ) 
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•1.2  Apparent  Radiated  Power 

In  Figures  3  to  i)  the  dashed  lines  plot  I  F  •  x|  from  Kq.  (4)  for  a  flat  eartli  id 
conductivity  0.02  S/m  from  a  5 ti 0  kllz  transmitter  radiating  1  k\V.  Regression 
analysis  on  KAOC  measured  data  over  path  segments  where  no  mountain  ridges  ..re 
interposed  between  the  transmitter  and  observer  (for  example.  Figures  7  and  H) 
yields  y-intereepts  displaced  approximately  1.6  dll  below  the  equivalent  i  k\\  nominal 
\alue.  This  discrepancy  could  be  due  to  the  NJU-25T  calibration,  since  J.  ■  x  data 
supplied  by  K.HLL  is  consistently  higher  than  the  KADC  measurements.  On  the  other 
hand,  beyond  the  Gila  .Mountains  (Figures  3  through  6)  the  regression  lines  are  dis¬ 
placed  from  3  to  5  dll  below  the  nominal  1  kU  value.  This  reduction  in  apparent 
radiated  power  can  be  interpreted  as  the  effect  of  the  mountain  ridge  in  extracting 
energy  from  the  propagating  ground  wave. 

The  losses  may  be  due  to  scattering  by  the  terrain  elevation  changes  and/or 
to  increased  absorption  bv  a  poorly  conducting  path  segment,  but  neither  mechanism 
can  be  adequately  treated  by  flat  earth  theory.  Indeed  for  certain  paths  the  signal 
strength  drops  by  as  much  as  15  dll  behind  the  mountain  ridges  (Figures  6  and  i)  and 
then  exhibits  a  recovery  effect  with  distance.  Such  a  variation  requires  the  more 
comprehensive  theory  which  follows  in  the  next  section. 


5.  PROPAGATION  OVER  ROUGH,  NON-IIOMOGENEOES  TERRAIN 
5.1  lliilford's  Integral  Equation 

Calculations  of  the  effects  of  irregular  terrain  can  be  made  via  a  one-dinien- 

o 

sional  integral  equation  derived  by  llufford.  '  llis  Kq.  (1  1)  in  the  notation  of 
Figure  *)  is 


W(x>  - 


f  \V(s)  ‘  f(s)  -  ■■  ■■- 
o 


(7) 


where  U  is  the  attenuation  function  for  the  Hertz  potential.  In  terms  of  the  vertical 
electric  field,  an  approximate  expression  for  \\  is 


!).  llufford,  G.  A,  (1.052)  An  integral  equation  approach  to  the  problem  of  wave 
propagation  over  an  irregular  terrain,  Quart.  .1.  Appl.  Math.  !l;3!U-404. 
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provided  VV  is  slowly  varying.  The  quantity 


f(s)  =  e 


-  in  /4 


kT  \f  8r2  1  ik(rl+r2-r0) 

2T  1° +  5TTJ  ° 


where  5  is  the  (normalized)  surface  impedance 
-  i  7i  /  4 


{' 


Figure  it.  Geometry  lor  the  Integral  liquation, 
s  and  x  are  the  horizontal  distances  of  the 
scattering  point  Q  and  the  observer  P, 
respectively,  from  the  source  T.  ft  is.  the  slope 
of  r^,  p  is  the  terrain  slope  at  Q,  and  n  is  a 
unit  vector  at  Q  normal  to  the  terrain 


The  factor  dr^ldn  takes  into  account  the  terrain  slope  and  the  aspect  angle  with 
respect  to  the  observer.  As  may  be  seen  from  Figure  9, 

9r2 

=  sin  (p  -  O') 


where  0  is  the  terrain  slope  and  <v  is  the  slope  of  r... 

For  gently  undulating  terrain  with  gradual  inhomogoneitics,  the  functions  VV(s) 
and  f(s)  are  slowly  varying.  To  the  extent  that  VV  and  f  can  be  represented  as  con¬ 
stants  within  arbitrarily  small  intervals  s  =  s  -  s  an  approximate  solution 

n  n  -  t 

ot  Kq.  (7  )  is 


W 


n  + 


VV  f  I 
n  n  n 


(8) 


where  W  is  a  constant  value  for  VV  over  the  nth  interval,  f  is  a  constant  value 
n  n 

for  f  over  the  same  interval  and 
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5.2  Numerical  Solution  lor  Special  Cases 

Equation  (8)  was  solved  numerically  via  a  CDC  6600  computer  lor  several 
cases.  In  Figure  10  the  results  of  including  a  Gaussian-shaped  ridge 

G(x>  =  h  exp  £  -9  >2j  (9) 

on  an  otherwise  flat,  homogeneous  earth  are  shown  for  two  values  of  conductivity. 

In  Eq.  (9),  G(x)  is  the  terrain  elevation,  h  is  the  ridge  height  at  x  =  b  and  w  is  the 
ridge  width  measured  at  G(x)  -  h/10.  This  form  corresponds  to  the  model  chosen 
by  Berrv.  ^  The  attenuation  function  W|  at  first  decreases  at  the  flat  earth  rate, 
and  then  increases  to  a  maximum  just  before  the  crest  of  the  ridge.  Behind  the 
ridge  there  is  a  minimum,  and  then  a  partial  recovery  with  increasing  distance 
from  the  source. 

Figure  11  illustrates  the  perturbation  in  |w|  caused  by  a  6  km  wide  segment 
of  conductivitv  cr  ,  centered  at  80  km,  in  an  otherwise  homogeneous  flat  path  of 
conductivity  0.  02  S/m.  Such  a  path  represents  the  much  poorer  conductivity  ex¬ 
pected  in  the  mountain  ridges.  **  |w|  follows  the  flat  earth  attenuation  rate  out  to 

the  beginning  of  the  segment,  and  then,  depending  on  the  selected  value  for  ffg. 
decreases  more  or  less  rapidly  until  the  far  edge  of  the  segment  is  reached.  As  in 
the  case  of  the  ridge,  |  VV  |  partially  recovers  as  the  distance  from  the  inhomo¬ 
geneity  increases.  (The  abrupt  conductivity  change  in  this  model  does  not  satisfy 
the  conditions  under  which  Eq.  (8)  was  derived,  so  the  results  are  not  valid  inthe 
immediate  vicinity  of  the  edges  of  the  segment.  )  A  comparison  of  Figures  10  and  11 
with  7  and  8  seems  to  show  that  the  effect  of  a  segment  of  poor  conductivity  more 
closely  matches  the  observed  |e  -  x|  variation  than  the  elevation  effect.  The  models 
demonstrate  that  irregularities  in  both  profile  and  earth  conductivity  play  significant 
roles  in  determining  the  propagation  of  MF  over  ridges. 


10.  Berry,  L.  A.  (1967)  Radio  propagation  over  a  Gaussian-shaped  ridge,  IEEE 

Trans.  Antennas  Propag.  ^AP  -  15(No.  5):701-7  02. 

11.  FAO-Unesco,  Soil  Map  of  the  World,  Vol.  II,  Unesco-Paris,  1974. 
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GAUSSIAN  RIDGE  AT  >-32 
(h  -0.4km,  w  •  7  5km) 


Figure  10.  The  Attenuation 
Function  Wl  in  the  Surface 
>f  tlie  Midge  Shown  at  the 
Bottom  of  the  Figure,  for 
Two  Values  of  Karth 
Conductivity 
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Figure  11.  The  Attenuation 
Function  Iwl  >n  the  Surface 
of  the  Inhomogeneous  Path 
Shown  at  the  Bottom  of  the 
Figure,  for  Three  Values 
or  a 

s 


<T  -.02 


(>.  PK.NKTRATION  ItFLOW  TI1K  DKSI.RT  FLOOR 


Attenuation  of  All-'  fields  with  depth  i n  homogeneous,  well-conducting  soil  is 
described  by 


-v/6  -jv/5 

11-11  e  s  e  s 
o 

where?  11^  is  the  magnitude  of  the  field  at  the  surface,  y  is  the  depth  below  the  sur¬ 
face  and  6  .  is  the  depth  of  penetration  or  skin  depth.  Attenuation  measurements  of 
tne  magnetic  fields  penetrating  the  desert  soil  were  made  by  lowering  a  loop  an¬ 
tenna  into  18-in.  diameter  holes  drilled  approximately  6-ft  deep  bv  a  commercial 
power  auger,  and  recording  the  decibels  change-  from  the  surface  reading.  An 
effective  skin  depth  was  then  calculated  via 

6  -  8-69  y  n  m 

s  IT  |  dB]  -  TTJ  dB]  '  <10> 

A  total  of  .10  holes  were  drilled  at  various  locations  as  indicated  in  F  igure  1.  The 
relative  attenuation  was  read  to  within  about  ±  0.  25  dB  on  the  N1U-25T  panel  meter, 
and  the  depth  was  measured  accurately.  The  average  skin  depth  calculated  was 
28.  25  ft  with  a  standard  deviation  of  18.  00  ft. 


7.  DISCUSSION  AND  CONCLUSIONS 

The  data  in  this  report  can  be  used  to  estimate  total  propagation  path  loss  for 
a  point-to-point  communication  link  between  terminals  buried  in  the  desert  south¬ 
east  of  Yuma  and  separated  by  mountain  ridges.  The  desert  valley  soils  have  rela¬ 
tively  high  apparent  conductivity,  with  a  correspondingly  low  attenuation  of  propa¬ 
gating  MF  radio  waves.  The  mountain  ridges  can  have  a  large  effect  which  depends 
on  the  location  of  the  terminal.  The  mountains  can  be  modeled  to  predict  the 
attenuation,  but  good  accuracy  requires  that  both  the  terrain  profile  and  the  ground 
constants  be  known.  An  overburden  of  desert  soil  will  result  in  additional  loss. 

The  measurements  and  portions  of  the  modeling  work  in  this  report  were  pre¬ 
sented  to  an  MX  C  working  group  at  a  technical  interchange  meeting  at  Norton  AFB, 
California  on  14  December  1977. 
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